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ABSTRACT 

The intermittent behavior of the space-scale distribution of Lja transmitted 
flux of QSO HS 1700+64 has been analyzed via a discrete wavelet transform. 
We found that there are strong indications of intermittency on scales down to 
about 10 kpc. These are: 1.) the probability distribution function of the 
local fluctuations of the flux is significantly long-tailed on small scales, and 2.) 
the local power spectrum of the flux shows prominent spiky structures on small 
scales. Moreover, the local power spectrum averaged on regions with different 
sizes shows similar spiky structures. Therefore, the random mass density field 
traced by the Lya forests is rougher on smaller scales, consistent with singular 
clustering. 

Subject headings: cosmology: theory - large-scale structure of the universe 
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1. Introduction 

The structures of the cosmic mass field on scales from sub-Mpc to kpc have attracted 
a lot of attention recently. High resolution N-body simulations show that the core profiles 
of massive halos of the cold dark matter (CDM) cosmogony are singular (Navarro, Frenk 
& White 1997, Moore et al. 2000, Jing & Suto 2000), while the halo profiles required by 
the rotation curves of dwarf galaxies (Flores & Primack 1994, Burkeret 1995) are shallower 
than the numerical results. The central cusps of dark halos also disagree with soft halo 
profiles inferred from low surface brightness galaxies (de Blok & McGaugh 1997.) The 
cores of galaxies and clusters are even found to be consistent with the thermal equilibrium 
model with an "universal" mass density (Firmani et al. 2000.) In other words, the singular 
behavior of cosmic clustering has not been detected in the cores of galaxies and clusters. 

We can now ask the question: can the singular behavior, if it exists, be revealed by 
methods other than the mass profile of galaxies and clusters? At the first glance, this goal 
seems to be not achievable, as singular mass density profiles can only be seen in cores 
of galaxies and clusters. However, a random mass field p{x) consisting of rare singular 
structures randomly scattered in low mass density background typically is intermittent 
(Zel'dovich et al. 1990.) A basic character of an intermittent field is that the density 
difference between two neighboring positions, \p{x + r) — p{x)\, can be "abnormally" large 
when r is very small. That is, the rare events of large density difference \p{x + r) — p{x)\ 
on small scales r have a higher probability than that for a Gaussian field. Such a singular 
behavior of a random mass field means that the probability density function (PDF) of the 
density differences on small scales, r, has a long tail. Obviously, the effects of long tailed 
PDF of the density difference are not limited to the singular mass profile. This motivated 
us to look for the PDF's long tail and its effects by using samples other than the cores of 
galaxies and clusters. 
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The PDF's long tail have not yet been seriously studied. The most popular statistical 
measure of large-scale structures - the power spectrum of a mass field, is insensitive to 
the PDF's long tail. Furthermore, The density difference, \p{x + r) — p{x)\, is a quantity 
localized in space x and on scale r, and so a space-scale decomposition is necessary. 
Thus using the power spectrum or any statistic that is not based on proper space-scale 
decomposition, it is not possible to identify the effects of long-tailed PDF of the density 
difference. 

In this Letter, using a discrete wavelet transform (DWT), we look for the long tail 
effects from a sample of QSO's Lya forests. It is believed that the distribution of baryonic 
diffuse matter is almost point-by-point proportional to the underlying dark matter density. 
Moreover, the absorption optical depth of Lya is linearly dependent on the baryonic density. 
Therefore, the high resolution data of the transmitted flux of QSO's absorption would be a 
good candidate to reveal the long-tailed PDF of cosmic mass field on small scales. 



2. Method 

Let us consider a 1-D random mass density field p{x) in spatial range L. With a DWT 
space-scale decomposition, the local density difference, \p{x + r) — p{x)\, is represented by 
the wavelet function coefficients (WFCs) as 

~^3,l = {■^jhP)^ (1) 

where ipj^i{x) is the orthonormal and complete basis of the discrete wavelet transform, and 
(...) = / ....dx is the inner product (Daubechies, 1992.) We use DAUB4 wavelet (Press et al. 
1993, Nielsen 1998) for our analysis through out this Letter. The WFC, e^y, is the density 
fluctuation on the scale L/2^ at the position i = 0, ...2-' — 1, or the mean density difference 
between nearest neighbors ranging on scale L/2^ at /. If the "fair sample hypothesis" 
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(Peebles 1980) holds, then the 1^ values of Ij^i form an ensemble of the density differences 
on scale j, and therefore, the distribution of e^^ is a reasonable estimate of the PDF of the 
density differences on scale j (Fang & Thews 1998.) 

The second order statistics je^yp describes the power of the perturbations of the mode 
(j, /). In other words, at a given position /, the local power spectrum is given by 




By averaging P^^x over all positions /, we have 

2^-1 2^-1 

It has been shown that Pj actually is a band-averaged Fourier power spectrum (Pando & 
Fang 1998; Fang & Feng 2000.) 

The Fourier power spectrum lacks phase information, and therefore, Pj cannot show 
the phase-related features of clustering. However, P^y is phase-sensitive. One can search for 
the phase-related features of the mass field by local DWT power spectrum P^y. 

We can generalize the definition of local DWT power spectrum, eq.(2) as follows. First 
we chop L into 2^' sub-interval, labeled by Is — 0, l...(2^^ — 1). Each sub-interval has a 
length LjT" . Then, the local DWT power spectrum at sub- interval Ig is given by 

Pj,{j.,is} = ^pj: 12 I^jvT- (4) 

l=ls2i-is 

It is the power on scale localized on Ig with size L/2■^^ 

For a Gaussian field, the local power spectrum Pj i will not show structures with 
respect to /. On the other hand, the singular behavior of a random field is measured by the 
exponent a defined by \p{x + r) — p{x) \ ~ r". The larger the a is, the smoother the field on 
small scales is, and vice versa. If the exponent a is negative there are an actual singularity 
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of the field. Therefore, the singular behavior can be revealed by the roughness of the local 
power spectrum on small scales. The WFC local power spectrum can also measure the 
index n of power-law profile p ~ r~" for individual core. 

3. Sample and analysis 

The sample used for the analysis is the Lya transmitted fiux of QSO HS1700-I-64. 
This sample has been employed to study the evolution of structure (Bi & Davidsen 1997), 
the Fourier and DWT power spectra (Feng & Fang 2000.) The recovered power spectrum 
has been found to be consistent with the CDM model on scales larger than about 0.1 
Mpc. The data ranges from 3727.012A to 5523. 554A with a resolution of 3 kms~^, for 
a total of 55882 pixels. In this paper, we use the first 25000 pixels for analysis, which 
correspond to ^ = 2.07~ 2.65, or A = 3727.012A ~ 4434.266A. On average, a pixel is 
about O.O29A, equivalent to physical size ~ 5 kpc at 2; ~ 2 for a fiat universe. We pad 
7768 null pixels at the end to utilize a fast wavelet transform algorithm, which requires 
the data size in powers of 2. It does not affect the analysis, because the wavelet transform 
is localized. Moreover, we subject DWT directly to pixels without transforming them to 
physical positions. The relation between pixel number and physical position is not linear, 
but it does not affect structures on small scales. Thus, we ignore the effect of the non-hnear 
relation in our present analysis. 

Most lines in the Lya transmitted fiux are due to absorptions by gases in cool and 
low density regions. The pressure gradients are generally less than gravitational forces. 
That is, the gas, and hence the transmitted fiux, should be good tracer of the dark matter. 
Nevertheless, small scale structures of the dark matter field may be smoothed out by the 
velocity dispersion of Lya forest gases. Therefore, to identify the clustering feature, we 
will statistically compare the real data with its phase-randomized counterpart, which is 
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obtained by taking the inverse transform of the Fourier coefficients of the original data after 
randomizing their phases uniformly over [0, 2n] without changing their amplitudes. 

3.1. The PDFs of WFCs 

In Fig. 1, we show the PDFs of the WFCs for j = 8, and 14. Each PDF is normalized 
to have unit variance. For scale j — 8, the departure from Gaussian distribution is not so 
significant. Especially, no tail shows in the j — S PDF, i.e. no WFCs found to be ij^i > 3a. 
For j = 14, the PDF of WFCs has 2^^ events. Therefore, if the PDF were Gaussian, the 
number of the events larger than 3a would be about 44. However, the data shows 234 
events beyond the 3a range. Furthermore, a Gaussian PDF predicts that the number of 
the events larger than 5a should be 0.01, while the data shows more than 100 events larger 
than 5a. The data extends to beyond 15a on both sides. Therefore, the PDF is indeed 
significantly long-tailed on small scales. In other words, the field is rougher on smaller 
scales. This indicates that the field may contain singular structures. 

The shape of the two PDFs of j = 8 and 14 is very different from each other. That is, 
the two stochastic variables e^,;, j — 8, 14 don't relate to each other as 

6,, = 2^(^-^')6,,„ (5) 

where /3 is a constant. Therefore, the mass field traced by the QSO HS1700-I-64 is unlikely 
to be self-similar. 

3.2. Local power spectrum 

In Fig. 2 we plot the local DWT power spectra of the HS1700-I-64 transmitted fiux and 
its phase-randomized counterpart. We take jg = 7, i.e. chopping the entire sample into 128 
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sub-intervals, and, in each sub-interval, calculating the power spectra for j' = 8 ~ 14, which 
correspond to physical scales 2^^~^ x 5 kpc. 

Fig. 2 shows that the j — S local power spectrum for real data is not very different 
from its phase-randomized counterpart. This is consistent with the j = S PDF shown 
in Fig.l. It is closer-to-Gaussian. While the j = 12 ~ 14 local spectra are very rough, 
showing remarkably spiky structures, that completely disappear in the phase randomized 
counterpart. The spiky features mean that a significant part of the power is concentrated in 
some small areas. This feature is a result of the long tailed PDF of the density difference, 
i.e., higher probability of "abnormal" density change. The smaller the scale, the more 
pronounced the spiky features. This, again, points to a singular clustering of cosmic mass 
field. 

It should be emphasized that the spikes in the local power spectrum with high j do 
not always correspond to the peaks in the density distribution (or the absorption lines in 
the optical spectrum). The WFCs e^y describe the difference in density between intervals of 
length 2^^~^ x 5 Mpc. The average of e^y over I generally is zero. The mean power (or 
variance) at j — 14 is Pu — 1.2 x 10~^. Thus, even a single event > 10a at j — 14 doesn't 
always refer to high density, and it can happen in regions other than high density cores. The 
spikes denote the positions where the density (or the absorption optical depth) undergoes a 
dramatic change, which is the key indicator of the singular behavior of a random field. 

The wavelet functions -0^^; are orthogonal, and therefore, the local power spectrum on 
scale j doesn't contaminate perturbations on other scales. This is very different from the 
density distribution smoothed by a window function on scale j. The peaks identified from 
the window-function smoothed density field contain all contributions from perturbations on 
scales / < j. Therefore, the peaks in a smoothed field are actually given by a superposition 
of perturbations on large and small scales. They may not show singular features, because 
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the PDFs of large scale (or / < j) perturbations are closer to Gaussian. 

We should estimate the possible distortion of the long tail effects caused by velocity 
dispersion of gases. We calculated the j — 14 local power spectrum with sub-interval 
js = 14, which is plotted in Fig. 3. This local power spectrum has almost the same spiky 
features identified in the j = 14, jg = 7 spectrum (see. Fig. 2.) That is, most spikes 
shown by the power localized in sub-interval with size of about 600 kpc actually are 
localized in sub-interval with size only about few 10 kpc. This result indicates that 
the contamination of gas velocity dispersion may not be significant, at least, for prominent 
spikes. Thus, the spiky structures should mainly come from the underlying mass field. 

Intermittcncy can more clearly be seen in Fig. 3. The mean powers [cq.(3)] of the real 
data (left panel) and randomized counterpart (right panel) of the js — 14 and j — 14 local 
spectrum actually are the same, while the spikes of the real data are higher than the mean 
power by a factor of few tens even hundreds. That is, in the real sample, most power of the 
j — 14 perturbations is concentrated in the spikes, and almost no power, i.e. Pj,{j^,i^} — 0, 
in places other than the spikes. This is a typical intermittent distribution. 

4. Conclusions 

With the PDFs of ij^i and local DWT power spectra, we show that the mass field traced 
by the QSO HS1700-I-64 Lya forests is neither Gaussian, nor self-similar, but intermittent. 
The spiky features shown in the local DWT power spectrum is remarkably pronounced 
on scales down to about 10 kpc. Moreover, the long tail and the spiky features are 
substantial on smaller scales. This indicates that the cosmic mass field is rougher on smaller 
scales, consistent with singular clustering. 

A big advantage of intermittcncy is that one can detect singular clustering using the 
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st atistical features of entire random density field, not limited to the cores of galaxies and 
clusters. The information of intermittency extracted from Lya forests would be important 
to test models of cosmic clustering in terms of their singular behavior. 

We thank Dr. D. Tytler for kindly providing the data of the Keck spectrum HS1700+64. 
PJ would also like to thank Dr. Robert Maier for his help. HZ thanks Dr. David Burstein 
for helpful discussions. 
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Fig. 1. — The PDFs of the WFCs e^y for j — 8 and j — 14 are shown on top right corner 
of each panel. The tail is shown on the main part each panels. Horizontal axis is for ej^i/a, 
where a is the variance of the sample. Vertical axis is the probability density. The solid 
curve is the Gaussian distribution with zero mean and unit variance. The PDF of j — 14 at 
zero is ~4.5, so it is not shown in the figures. 
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Fig. 2. — The vertical axis represents '^^~^''Pj,{js,is} with jg — 7. The left and right panels 
are for the HS 1700+64 transmitted flux, and its phase-randomized counterpart, respectively. 

The horizontal axis represents the position Ig in units of 2^ x 10 h^^ kpc. The scales of the 
local power spectra, j = 8, 10, 12, 14, correspond to physical scales 2^^~^ x 5 kpc. 
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Fig. 3. — The same as Fig. 2 with j = 14, but taking jg = 14. 



